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Momentum Dependent Local-Ansatz Approach to Correlated Electron 

Systems: Non Half-Filled Case 
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Momentum dependent local-ansatz wavefunction approach (MLA) to the correlated elec- 
tron systems in solids has been further developed to solve best a self-consistent equation 
for variational parameters at non half-filling. With use of the improved variational scheme 
we performed the numerical calculations for the non-half-filled band Hubbard model on the 
hypercubic lattice in infinite dimensions. We verified that the self-consistent scheme signif- 
icantly improves the correlation energy and the momentum distribution as compared with 
the original scheme in the MLA. We also demonstrate that the theory improves the standard 
variational methods such as the Local-Ansatz approach (LA) and the Gutzwiller wavefunc- 
tion approach (GA); the ground-state energy in the MLA is lower than those of the LA and 
the GA in the weak and intermediate Coulomb interaction regimes. The double occupation 
number is shown to be suppressed as compared with the LA. Calculated momentum dis- 
tribution functions show a distinct momentum dependence, which is qualitatively different 
from those of the LA and the GA. 

KEYWORDS: variational method, electron correlations, Gutzwiller wavefunction, local ansatz, 
Hubbard model, critical Coulomb interaction, infinite dimensions 



1. Introduction 

Electron correlations play an important role for understanding the electronic structure, 
metal-insulator transition, and the high-temperature superconductivity in the solid-state 
physics. Thus many theories have been proposed so far to describe correlated electron sys- 
tem 1,2 ) on the basis of the variational method, 3 ~ 5 ) the Green function techniques, as well as 
many numerical techniques such as the Monte-Carlo method. 6 ~ 8 ) 

The variational theory among various methods has been developed as a practical method 
for understanding the ground-state properties of correlated electrons from molecules to solids 
over 50 years. A minimum basis set to describe correlated electrons is constructed in this 
approach by applying one-particle, two-particle, and higher-order particle operators onto the 
Hartree-Fock wavefunction, and their amplitudes are chosen to be best on the basis of the 
variational principle. The Gutzwiller wavefunction is one of the popular trial wavefunction 
in solids because of its simple and intuitive structure. This approach was first introduced by 
Gutzwiller to clarify the role of electron correlations in metallic ferromagnetism. 9-11 - 1 There, 
one reduces the amplitudes of doubly occupied states on the local orbitals in the Hartree-Fock 
wavefunction by making use of a projection operator IT(1 — gn^nii). Here nj CT is the number 
operator for electrons on site i with spin a. The variational parameter g is determined by 
minimization of the ground-state energy. Brinkman and Rice recognized that the Gutzwiller 
approximation describes the metal-insulator transition. 12 ) Because the Gutzwiller method is a 
nonperturbative approach, it has extensively been applied to the strongly correlated electron 
systems. 13 ^ 

The Gutzwiller wavefunction in the Gutzwiller ansatz approach (GA) yields a physical 
picture of electron correlations and is useful for correlation problems, but it was not so easy 
to apply the method to realistic Hamiltonians. The approach was successfully generalized 
by Stollhoff and Fulde 14 ~ 16 ) by using an alternative method called the local-ansatz approach 
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(LA), which is simpler in treatment and applicable to realistic Hamiltonians. The LA takes into 
account the excited states created by local two-particle operators such as {Oi} = {Sn^Snn}, 
and determines their amplitudes variationally. Here 5rii a = rii a — (nj CT )o, (^mt)o being the aver- 
age electron number on site i with spin a in the Hartree-Fock approximation. The theory has 
been applied to many systems such as molecules, transition metals, polyacetylene, transition 
metal oxides and semiconductors. 1, 17 - ) 

Although the LA is able to explain fruitfully the correlation effects in actual materials, 
the application has been limited to the weakly correlated region because of the difficulty in 
evaluation of the higher-order terms in average quantities. The expansion of the Hilbert space 
by the local operators is not sufficient to characterize precisely the weakly correlated states; 
the LA does not reduce to the second-order perturbation theory in the weak correlation limit. 
For the Gutzwiller wavefunction, the same difficulty also arises even in infinite dimensions. To 
overcome the difficulty, Kakehashi et. a/. 18 ) proposed a variational wavefunction theory called 
the momentum-dependent local ansatz approach (ML A). 

When we expand the local operators {Oi} in the LA by means of the two particle oper- 
ators in the momentum space, we find that each coefficient of the expansion is momentum 
independent. In the MLA wavefunction, 1S ) we replace the constant coefficients in the LA with 
the momentum-dependent variational parameters in order to obtain the best local operators. 
It results in a new set of local operators {Oi}. We then construct the MLA wavefunction with 
use of the local operators {Oi} as I^mla) = — Oijlfo)- Here |</>o) is the Hartree-Fock 

wavefunction and i denotes site of atoms. The best local basis set is chosen by controlling 
the variational parameters in the momentum space. We calculate the ground-state energy us- 
ing the MLA wavefunction within a single-site approximation (SSA). Minimizing the energy, 
we obtain a self-consistent equation with variational parameters. It is however difficult to 
solve the self-consistent equation directly. Because of this, we obtained in the previous paper, 
which we refer to I, an approximate solution which interpolates between the weak Coulomb 
interaction limit and the atomic limit. 

In this paper we point out that it is indispensable toward quantitative calculations to 
choose the variational parameters best, though in our paper I 18 ) we applied approximate vari- 
ational parameters, and improve the variational parameters on the basis of the variational 
principle. We investigates the validity of our theory for the non half-filled case performing 
numerical calculations of various physical quantities. Especially for non half-filled case we ob- 
serve that the best choice of variational parameters gives reasonable results for the momentum 
distribution, while the previous version of the variational parameters yields unphysical results 
near the Fermi level. Moreover, we demonstrate that the improved variational parameters 
much improve the LA in the weak and intermediate correlation regimes. 

The outline of the paper is as follows. In the following section we write down our wavefunc- 
tion for the single-band Hubbard model. We obtain the ground-state energy within the SSA 
and derive the self-consistent equation for the momentum dependent variational parameters 
on the basis of the variational principle. We develop the theory to obtain the best value of 
variational parameters to solve the self-consistent equation. In §3, we present our results of 
numerical calculations. We will clarify the role of the best choice of variational parameters 
on various quantities. Furthermore, we discuss the correlation energy, the double occupation 
number, the momentum distribution function, and the quasiparticle weight as a function of 
the Coulomb interaction energy parameter, and verify that the present approach improves the 
LA in the weak and intermediate Coulomb interaction regimes for non half-filled band. We 
summarize our results in the last section and discuss the remaining problems. 

2. Momentum-Dependent Local Ansatz with the best Variational Parameters 

We adopt in this paper the single-band Hubbard model 18-22 ) as follows. 



H = ^2(e - <rh)n ia + ^ Uj a\ 




(1) 
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Here eo (h) is the atomic level (magnetic field), ty is the transfer integral between sites i and j. 
U is the intra-atomic Coulomb energy parameter. a\ a (aj CT ) denotes the creation (annihilation) 
operator for an electron on site i with spin a, and = a\ a ai a is the electron density operator 
on site i for spin a. 

In the Hartree-Fock approximation, we neglect the fluctuation term and replace the many- 
body Hamiltonian (1) with an effective Hamiltonian for independent particle system, 

H o = Yl ^ a ^ a i° ~ U Y1 ( n it)o{na)o , (2) 

ija i 

and approximate the ground-state wavefunction |^) with that of the Hartree-Fock Hamil- 
tonian Hq, i.e., \4>o). Here Uj a = (eo + U{rii- a )o — crh)5ij + Uj(l — 5ij). (~)o denotes the 
Hartree-Fock average (<^o|(~)|0o)> and (nj (T )o is the average electron number on site i with 
spin a. The original Hamiltonian (1) is then expressed as a sum of the Hartree-Fock Hamil- 
tonian (2) and the residual interactions as follows. 

H = H + uY,O l . (3) 

i 

Here Oi = Sn^Sn iir and 5n ia = n ia - (n ia ) . 

In the LA, 16 ) the Hilbert space expanded by the local operators such as the residual 
Coulomb interactions {Oi} = Jn^Jn^ is taken into account in order to describe the weak 
Coulomb interaction regime. The ansatz for the Hubbard model is written as 

l*LA>= [II^-^aO*)]!^ • (4) 

i 

Here r/ LA is the variational parameter as the amplitudes of the basis set expanded by {Oi}. 
The LA is different from the Gutzwiller ansatz wavefunction |^ga) = Yli(^ ~ 

gni^rii^) \<p Q ) in which the doubly occupied states are explicitly controlled by a variational 

parameter g (0 < g < 1), and simplify the evaluation of the physical quantities in the weakly 
correlated region. 

As we have emphasized in our previous paper I, 18 ) though the LA is useful for under- 
standing correlation effects in real system, the Hilbert space expanded by the local operators 
{Oi} is not sufficient to characterize exactly the weakly correlated region; it does not reduce 
to the second-order perturbation theory. In order to describe the weak Coulomb interaction 
regime correctly we introduced a new set of local operator 

Oi= <*a l*> <*l fc i> <*2 1*> <*l*2>^fc a *ifc 1 <y(a^ Jr afc 2 4.)*(ai/ t Ofcit) i ( 5 ) 

kik2k[ k' 2 

and proposed the following new wavefunction with momentum dependent variational param- 
eters {r} k ' 2k2k ' lkl }. 

i*> = na - <wo> ■ (6) 

i 

Here (i\k) = exp(— ik ■ Ri)/y/N is an overlap integral between the localized orbital and the 
Bloch state with momentum k, Ri denotes the atomic position, and N is the number of sites. 
a \ a (flfco-) denotes the creation (annihilation) operator for an electron with momentum k and 

spin a, and 5(a\ la a k(T ) = a\> a aka ~ (a\> a ako)o- 

The operator Oi is still localized on site i because of the projection {k[\i) {i\ki) {k' 2 \i) (i\k2) ■ 
It should be noted that 0\ / Oi and OiOj / OjOi (i / j) in general. The wavefunction 
\ty) reduces to |^la) when {rj^'^k^ki} become momentum-independent because Oi — > t/^aOi 
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when %^ 2fc / fcl -»• rj LA . 

The variational parameters are determined by minimizing the ground-state correlation 
energy E c . 

E c = (H) - (H) = i±±2 . ( 7 ) 

Here H = H — (H) . 

Although it is not easy to calculate the correlation energy with use of the wavefunction 
(6), one can obtain the energy within the single-site approximation (SSA). In the SSA, the 
average of (A) of an operator A = A — (A)q with respect to the wavefunction (6) is given as 
follows: 

i ((i-ojxi-cy),, 

The derivation of the above formula has been given in Appendix A of our paper I. 18 ) 

By making use of the above formula, one can obtain the correlation energy per atom as 
follows. 

-(OlH)o - (Hdi) + {d\Hdi) 



l + (0?Oi)o 



(9) 



Here we assumed that all the sites are equivalent to each other for simplicity and we made 
use of the fact (Oj) = (6j)o = 0. 

Each term in the correlation energy (9) can be calculated by making use of Wick's theorem 
as follows. 

(Hdi) = u Y, <*;il*><*|fei)<^|i><*|fe>^<^|j)01^i><feU)0'l^2) 

x^fc 2 fc' 1 fc 1 /(4 1 t)( 1 - /(4it))/(^24)( 1 - f(h' 2 i)) > (10) 
(6}H)o = (Hdi)o » (11) 

(o}Hd t )o = £ m)(k^)m)(k 2 \i) V i /2k2k , ki 

k\k2k' 1 k' 2 

x/(e fclt )(l - /(eV lt ))/(e fe2 +)(l - /(e^)) £ (4lO(^ 3 )(^K)(^ 4 > 

^3 ^4 ^4 



x \^ E k l 2 k2k' 1 k^k 1 h.ih' 1 k'^k 2 k A h' 2 k' A + Uy^y^y^y^jriy^y^ , (12) 

U k ' 2k2k > lklk > 4kik ' 3k3 = uY\[{j\k 1 ){k z \j)f{l k ^)5 k > i y 3 - (fcl|j)0'l*3>(l - f(.h' 3 t)) S kik 3 ] 

j 

x[{j\k2){h\j)f(i k4i )5 k ' 2k ' 4 - (k' 2 \j){j\k' 4 )(l - f(e K x))Sk2kA , (13) 
(OjOi)o = ^4 E l^ 2 fcifcilVfet)(l-/(4 it ))/(^)(l-/(^|)) • (14) 

Here e fe(T = e ka - fi, e ka = e + U (n;_ CT ) + e k - ah and AE^^ = e k ^ - e k2i + e fe / t - e felt 
is a two-particle excitation energy. e& is the the Fourier transform of tij. 

The above expressions (10) and (13) contain nonlocal terms via summation over j (i.e., 
Ylj)- We thus make additional SSA called the R = approximation. 23 ' 24 ) For example, we 
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have in eq. (10) 

2^i|t><t|fci><^l<><<l*!i><*i|j>0"l*4><fe|j>0"l^> = l^e^-^W-*) , (15) 

3 j 

but we only take into account the local term (j = i). In the R = approximation, (HOi)o(= 
{d\H)l), and (d\HOi) Q reduce as follows. 

/(^it)( 1 -/(^t))/( g fe24)( 1 -/(efc^))%ife 2 fe' 1 fe 1 , (16) 



AT4 



+ ]V2 { /(^3t)/(^44-) Vk! 2 k 4 k[k 3 ~ ^ /(^3t)( 1 ~ f(^k' 4 l)) Vk'^kik^ks 



k 3 k 4 



X^ 1 ~ f ih'^)) f (h 4 \) Vk^kik'.^ + ~ /(^st))( 1 ~ /(%!.)) ^Ji^fci} 



Kg K 4 



(17) 



Variational parameters {f^fofe^} are obtained by minimizing the correlation energy e c , 
i.e., eq. (9) with eqs. (14), (16), and (17). The self-consistent equations for {"nk^k'^} in the 
SSA are given as follows. 

(Ai?jfc 2 jfc2fcifci - ^c)'']k! 2 k2k' 1 k 1 

U 



+ jyf [ f^k 3 -\)f{h 4 ,i)'nk l 2 k A k' 1 k3 - Y /(efc 3 f)( 1 ~ fih'^Vk'^k^ 

k 3 k' 4 



k 3 



k 3 k4 



- X^ 1 ~ f^k'^f^k^Vk'^k! + ~ /( g fc 3 t))( 1 - /(4^))^fc 2 fc 3 fciJ = ^ • ( 18 ) 

It should be noted that the variational parameters {% 2 'fc 2 fc 1 'fci} i n ec L- (18) vanish when 
U — > 0, i.e., r/k^k^ki ~ 0(U). Thus in the weak C7 limit, one can omit the second term at 
the l.h.s. (left-hand-side). We then obtain the solution in the weak U limit as 

U 

Vk'k2k',ki = it; • (19) 

In the atomic limit the transfer integrals tij disappear, and one electron energy eigen 
value efc becomes A;- independent, eo- Thus, AE k i k2k i kl vanishes. In this limit we can drop the 
k dependence of Vk'^k'^ii i- e -> r ?fe 2 fe 2 fejfci — ^ V- Then, we find a /c-independent solution being 
identical with the LA. 



■(OiHOi) + y/iOiHOift + ^OiHftiOfy 



" LA " 2(O t H) (Of) • (20) 

It is not easy to find the solution of eq. (18) for the intermediate strength of Coulomb in- 
teraction U. We therefore proposed in our paper I an approximate solution which interpolates 
between the weak and the atomic limits; we approximate {rj k > k 2 k' 1 k 1 } m the second term with 
the momentum-independent parameter 77 which is suitable for the atomic region. We have 
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(21) 



In the previous paper I 18 ) we made use of that in the LA for 77, and adopted the correlation 
energy in the LA for e c in the denominator. We call this the non-self-consistent MLA in the 
followings. 

The best value of rj, however, should be determined variationally in general. In this paper 
we further develop the theory in which 77 is determined best. According to the variational 
principle, the ground-state energy Eq satisfies the the following inequality. 

(mm 



E < E[V] 



(22) 



Here ^ is a trial wavefunction. 

In the MLA, we choose the wavefunction \I> = <5 {{Vk^k'^i}) an d the corresponding energy 
EdVk'^k'^}) satisfies the inequality E < E^rj^y^}) . For the stationary values vl'^k'^ 
, we have 



E < E{{r,* KkaKkl }) < E({ Vk , k2k[kl }) . 
In the previous calculations, 18 ) we obtained an approximate Vk'k 2 k' kl i^)> 

Ufj 



r )k i 2 k 2 k> 1 k 1 {'ni e c 



^E k ^ k2 y ikl - e c 



(23) 



(24) 



which is determined by the correlation energy e c and a momentum independent parameter fj. 
Here 

fj = [1 - r?(l - 2<n it )o)(l - 2(n 4 ) )]. (25) 
When we adopt the form (24) as a trial set of amplitudes, we have an inequality as 

(E ) < E({ V * k , k2k , ki }) < E({r) k , k2k , kl (fj,e c )}) . (26) 

The above relation implies that the best fj is again determined from the stationary condition 
of the trial energy E({r/ k2k2k ^ kl (fj,e c )}). Because e c should satisfy the stationary condition 
5e c = for the value fj*, fj* is determined by the following condition 

de({rj k2k2k/ikl (fj,e c )}) 



dfj 



= . 



(27) 



The self-consistent equation is obtained from eq.(27) in the same way as in eq. (18) 
7^4 E /( g fcit)[ 1 - /( g fc , 1 t)]/( g fc2l)[ 1 - f{h' 2 i)} ^ g ~ 1 1 (A-Efc/fc 2 fc^i 



k\k' x k2 k' 2 



k' 2 k2k' 1 ki 



fa\ ^2 ^2 



/(efe 3 t)/(^44-)%ifc4fe^ 3 ~ E^ ~ /(^t)]/(^44-) 7 ?^fc 4 ^fcl 



kgk4 



k'„k4 



E /(^3t)t 1 - fih'^Vk'^kiks + Et 1 ~ /( g fcst)][ 1 ~ f(h' 4 i)}vk' 4 k 2k ' 3 ki 

k 3 k' 4 



k' 3 k' 4 
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dilk'^k'^! U 
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Substituting the above expression into the self-consistent equation (28) we obtain 

1 

V = 



1 + 



uc_ ■ 
IT 



Here 



/(e fclt )[l - f(e k > lt )]f(e k2l )[l - /(e^)] 



kik' 1 k2k' 2 



(^^fe 2 ik' 1 fci ~ e c) 



{E 



E 

k',,k4 



[1 ~ /gfc|t)]/(ffc4jj 

{^Ek'kik'ki - e c) 



^ (A^fc/^fc/^ - e c ) ^ (A% 2 ^-e c ) J' 



and 



1 

AT4 



E 



/(e fclt )[l - /(e fcit )]/(e^)[l - /(e^)] 



- e r 



(A-E'^fc 2 /c / 1 /c 1 

In the energy representation, each term is expressed as follows: 

/(e~i t )[l - /(e2 t )]/(e 3 i)[l - f(hi)]f{hdf(hi) 



= 1 
-I 
-I- 



6 

11 de n 

n=l 



6 

EI P( e r, 
n=l 



6 

11 de r, 
n=l 



(e 4 - e 3 + e 2 - ei - e c )(e 4 - e 6 + e 2 - e 5 - e c ) 

/(?l t )[l " /(€2t)]/(?3l)[l " /(e4l)][l - /(e>)]/(^) 



n p(e 

n=l 



(e 4 — e 3 + e 2 — ei — e c )(e4 — e6 + 65 — e\ — e c 



6 

11 de r, 
n=l 



P(^n) 
n=l 



/(€»[! " /(e 2t )]/(63 4 )[l - /(€4l)]/(6 5t )[l - f(e 6l )] 



(e 4 - e 3 + e 2 - ei - e c )(e 6 - e 3 + e 2 - e 5 - e c ) 



+ 



and 



/ 



6 

EI de n 

n=l 



n P(e n ) 

n=l 



(29) 



(30) 



(31) 



(32) 



/(Slt)[l " /(€2t)]/(?3l)[l " /(€4l)][l " /(€5t)][l " /(^)] 



/ 



(e 4 - e 3 + e 2 - ei - e c )(e 6 - e 3 + e 5 - ei - e c ) 

f(e lt )[l-f(e 2t )]f(i u ) [1-/(641)] 



(33) 



4 
n=l 



P(*n) 
n=l 



• (34) 

e 4 — e 3 + e 2 — ei — e c 

Here e ncr = e n + e a and = e + U (n_ CT )o — ll is the atomic level measured from the chemical 
potential, and p(e) is the density of states for the one-electron energy eigen values for the 
non- interacting system i^. 
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It should be noted that the self-consistent solution (30) is also obtained by solving approx- 
imately the original self-consistent eq.(18). In order to do this, first we divide the both sides 
of eq.(18) by {/S.E k i^ k2 y^ kl — e c ) , substitute the form (24), and we obtain eq. (30) after taking 

the average with respect to hk'^k^ with a weight /(e felt )[l - f{e k '^)]f{h 2 i)[ l ~ f{h' 2 l)]- 
The variational principles on rj tells us that such a solution should be the best among possible 
approximate solutions. We also note that an approximate form (25), which was obtained in 
paper I, is derived from the solution (30). In fact, we rewrite eq. (30) as fj = 1 — UfjC/D. By re- 
placing the approximate form (24) in the expression of UfjC with the momentum independent 
value rj, we reach eq. (25). 

The ground-state correlation energy is obtained by substituting the variational parameters 
(24) into eq. (9). The each element in the energy is given as follows. 

(Hdi) = {d\H)i 



(6|#(A>o 
-u>?j 



4 

EI de n 

n=l 



Ln=l 



/(eit)[l-/(e 2 t)]/(€3i) [l-/(e4i) 



£4 — 63 + e 2 — £1 — e c 
(d\Hdi) Q = (OjfloOi)o + UidjOiOJo , 



4 

n de r. 
n=l 



4 

II P( e r, 
n=l 



f(ht)[l-f(^)]f(hi)[l-f(hi)} 



(e 4 - e 3 + e 2 - ei - e c ) 2 (e 4 - e 3 + e 2 - ei) 1 



u- 



4 

EI de r . 

n=l 



P(*n) 
n=l 



/(6 lt )[l-/(e 2t )]/(63 4 )[l-/(e 4 ;) 



(e 4 - e 3 + e 2 - ei - e c 



/ 



6 

. n=5 



ft P(tn) 
_ n=5 



/(e5t)/(eei) 



-/ 

-/ 



6 

. n=5 



(e 4 — €6 + e 2 — £5 — e c ) 



6 

n=5 



6 

EI de n 

. n=5 



(£6 - £3 + £2 - £5 - £c) 

[1 - /(£»]/(£6i) 



6 

11 P^r, 
n=5 



+ 



/ 



6 

11 de r . 

_ n=5 



(e 4 - £6 + £5 - £1 - £c) 

[1 "/(£»][! -fihi)] 



6 

EI P(£r, 
n=5 



(OT Oi ) = ^ 2 ^ 2 



/ 



A de r 

n=l 



(£6 - £3 + £5 - £l - £c) 
n P(£n)l /(£»[! - /(£»]/(£3i)[l - f{hi) 

n=l 



(35) 
(36) 



(37) 



(£4 - £3 + £2 - £1 - £ c 



(38) 



(39) 



It should be noted that fj in eq. (30) is given as a function of e c , and e c in eq. (9) depends 
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on fj and e c . Therefore, both equations have to be solved self-consistently. To determine the 
best value of fj, we start from e c in the LA for example, and calculate fj according to eq. (30). 
Next we calculate various elements (HO^q {{0\h)q), (OjHOi}^, and ((5jOj}o which are given 
by eqs. (35), (36), and (39), respectively. Using these values we calculate e c according to eq. 
(9). We repeat this cycle until the self-consistency of e c and fj is satisfied. We call this scheme 
the self-consistent MLA . 

Electron number (nj)(= Y2a( n iv))i the momentum distribution {n ka ), and the double 
occupation number (n^n^) are obtained from d(H)/deo, d{H)/de ka , and d(H)/dUi, respec- 
tively. Here e k a = £fc — ch. Making use of the single-site energy (9) and the Feynman-Hellmann 
theorem, 25 ) we obtain the following expressions. 



(rn) = (n;) + 
(nka) = (nka)o + 



it, 



{n^n iSr ) = (n it )o(n 4 ) + 



l + (OlOi) 

-{d\Oi) - (o t di)o + {o}oA)o + EAm-Mdlniaddo 



l + <oJOi)o 
Here hi = m - (nj) , h ka = n ka - {n ka ) , and 

(djfiMo 



5 

11 de n 

n=l 



P(<*) 
n=l 



f(h-a)[l ~ f{e 2 -a)]f{ha)[l ~ f (ha)] 



(e 4 - e 3 + e 2 - ei - e c ) 2 
[1 - f(ha)] f(hc 



(^5 ~ e 3 + e 2 - ei - e c ) (e 4 - e 5 + e 2 - ei - e c 



^<Otn fc<r Oi> 



[1-/(4.)] | 



3 

11 de r. 
n=l 



ft P(*n) 
n=l 



/(e 1 _ (7 )[l-/(e 2 _ CT )]/(e 3 ,) 



7(4.) | 



n 

n=l 



n 

n=l 



(efc<r - e 3 + £2 - ei - e c ) 2 

/(€i_ ff )[l ~ f(h-a)]f(ho 



(£3 - efca + e 2 - ei - e c ) s 



(40) 
(41) 
(42) 



(43) 



(44) 



(O]0i)o + (OiO^o 



= 2Ufj 



4 

n=l 



n P(e n ) 

n=l 



/(?it)[l-/(e 2t )]/(e 3 ^)[l-/(644)] 



(45) 



(e 4 - e 3 + e 2 - ei - e c ) 

Note that {OjOiO^o has been given by eq. (38). The expressions of these physical quan- 
tities consist of the multiple integrals up to the 6- folds. One can reduce these integrals up to 
the 2-folds using the Laplace transform. 26 ) Their expressions are given in Appendix. 

3. Numerical Results 

We have performed the numerical calculations for the non-half-filled as well as half-filled 
bands of the Hubbard model in order to examine the validity of the improved scheme of the 
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Fig. 1. The correlation energy E c vs. Coulomb interaction energy curve for n = 0.8. The thick solid 
curve: the MLA with self-consistent fj, the thin curve: the MLA with non self-consistent fj, and 
the dashed curve: the LA. 



MLA and the effect of electron correlations in the local ansatz. To calculate various physical 
quantities, we have adopted the hypercubic lattice in infinite dimensions, where the single-site 
approximation works best. The density of states (DOS) for non-interacting system is given 
by P( e ) = (Vv^) exp(— e 2 ) in this case. 18 ) The energy unit is chosen to be J dep(e)e 2 = 1/2. 
The external magnetic field h is assumed to be zero. 

3. 1 Role of the best choice of fj 

To calculate various quantities in the MLA, we solved the self-consistent equations (9) 
and (30) with use of the Laplace transforms of elements, which are given in Appendix. In this 
sub-section, we compare the self-consistent results with the non self-consistent ones to clarify 
the role of the best fj. 

Figure 1 shows the calculated correlation energy as a function of Coulomb interaction. The 
correlation energy for the MLA without best choice of fj gives the lower correlation energy as 
compared with the LA. The correlation energy for the MLA with the best choice of fj is lower 
than that of the non-self-consistent MLA. The results indicate that the self-consistency of fj 
is significant for finding the best energy. 

In Fig. 2 we show an example of the momentum distribution as a function of energy 
when electron number is deviated from 1. The MLA with non self-consistent fj (25) shows a 
bump in the vicinity of the Fermi level, leading to an unphysical result. The MLA with self- 
consistent fj yields a significant momentum dependence which shows monotonical decrease of 
the distribution with increasing e^ c . 

We have also calculated the quasiparticle weight Z vs. Coulomb interaction energy curves 
at half-filling. As shown in Fig. 3, we find that the best choice of fj increases Z (i.e., decreases 
the inverse effective mass), so that the critical Coulomb interaction of the divergence of the 
effective mass, U C 2 changes from 3.21 to 3.40. The latter is closer to the NRG 27 ) value U C 2 = 
4.10, which is believed to be the best at present. 

From the above discussions on the results with and without self-consistent 77, it is obvious 
that the best choice of fj improves the results of the MLA. In the following discussions we 
adopt the best choice of fj. 
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Fig. 2. The momentum distribution as a function of energy for various theories with n = 0.6 and 
U = 2.0. The solid curve: the MLA with the best choice of rj, the dashed curve: the MLA without 
the best choice of rj, and the dotted curve: the LA. 



1.2 




0.5 1 1.5 2 2.5 3 3.5 

U 

Fig. 3. Quasiparticle-weight Z vs. Coulomb interaction curves in the MLA with self-consistent r) 
(solid curve), and without (dashed curve). 



3.2 MLA vs LA in various physical quantities 

In this section, we present the numerical results on various physical quantities, and discuss 
the new aspects of the MLA and related electron correlation effects by comparing the MLA 
with the LA. 

In Fig. 4, we represent the calculated correlation energy per atom as a function of Coulomb 
interaction U . The energy in the MLA is lower than that of the LA over all Coulomb interaction 
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energy parameters U and electron numbers n. These results imply that the ML A improves 
the LA. The magnitude of the correlation energy |e c | tends to increase with increasing U, 
because with increasing U the correlation corrections increase as U 2 for small U and cancel 
the Hartree-Fock energy loss being linear in U for large U. For a fixed value of the Coulomb 
interaction U, the gain of the correlation energy |e c | increases with increasing n, because there 
is a correlation energy gain at each doubly-occupied site and the number of such sites increases 
with increasing n. 
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Figure 5 depicts the double occupation number (n^n^) vs. Coulomb interaction curves for 
the non-half-filled case. In the uncorrelated limit, the double occupancy is the same for both 
LA and MLA and it decreases with increasing Coulomb interaction U because electrons move 
on the lattice so as to suppress the loss of Coulomb energy due to double occupation. We find 
that the MLA wavefunction reduces the double occupancy as compared with that of the LA 
in the range < U < 5, while in the range 5 < U the double occupancy in the MLA is larger 
than that of the LA. It implies that the LA with momentum-independent tjla overestimates 
the itinerant character for weak and intermediate U regions, while it overestimates the atomic 
character for large U region. 

The momentum-distribution function shown in Fig. 6 indicates more distinct difference 
between the LA and the MLA. The distributions in the LA are constant below and above the 
Fermi level irrespective of U. The same behavior is also found in the GA. 9-11 ) The MLA curves 
show a monotonical decrease of the distribution with increasing efc CT , indicating a distinct 
momentum dependence of (n^) via energy e^, which is qualitatively different from both the 
LA and the GA. 




The quasiparticle weight Z (i.e. the inverse effective mass) is obtained from the jump 
at the Fermi level in the momentum distribution according to the Fermi liquid theory. 28 ' 29 ) 
Calculated quasiparticle weight vs Coulomb interaction curves are shown in Fig. 7 for the half- 
filled case. The quasiparticle weight in the LA changes as Z = (1 — 377 l A /16)/(1+ 77 la/16) an d 
vanishes at C/ c2 (LA) = 24/\/37r(= 7.82). In the GA, 12 ) the quasiparticle weight changes as 
Z = 1 — (U /U C 2) 2 ■ The curve in the GA agrees with the LA curve for small U. But it deviates 
from the LA when U becomes larger, and vanishes at U C 2(GA) = 8/y/ir (= 4.51). It should be 
noted that the GA curve strongly deviates from the curve in the NRG 27 ) which is considered 
to be the best. We observe that the critical Coulomb interaction f7 c2 for the self-consistent fj 
is 3.40 in the MLA while U C 2 in the non-self-consistent fj yields 3.21. The quasiparticle weight 
in the MLA much improves the LA as seen in Fig. 7. We note that the wavefunction itself 
does not show the metal-insulator transition at f7 c2 in the present approximation because 
the approximate expression of variational parameters (21) has no singularity at finite value 
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U 

Fig. 7. Quasiparticle-weight vs. Coulomb interaction curves in various theories. The RPT-1: dashed 
curve, 30 ) the NRG: thin solid curve, 27 ) the LA: dotted curve, the MLA: solid curve, and the GA: 
dot-dashed curve. 



of U. In this sense, the calculated Z and wavefunction are not self-consistent in the present 
approximation. The values of Z obtained by the LA and the MLA should be regarded as an 
estimate from the metallic side. 




Figure 8 also shows the quasiparticle weight as a function of Coulomb interaction energy U 
for the non-half-filled case and the half-filled case. In the uncorrelated limit, the quasiparticle 
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weight is 1 as it should be. It decreases with increasing Coulomb interaction U for both the 
MLA and the LA. In general, the curves in the MLA are lower than those in the LA as 
expected from the fact [/^(MLA) < U C 2(LA) at half- filling. This means that the electron 
effective mass is enhanced by the self-consistent treatment of fj irrespective of U and n. 

4. Summary and Discussions 

In the present paper, we improved the variational scheme of the MLA which self- 
consistently determines both the variational amplitude fj and the correlation energy e c making 
use of variational principles. To examine the improvement and validity of the theory, we per- 
formed the numerical calculations on the basis of the Hubbard model on the hypercubic lattice 
in infinite dimensions. We verified for both the half-filled and the non half-filled bands that 
the self-consistent scheme of the MLA improves the correlation energy, the momentum dis- 
tribution function as well as the quasiparticle weight. Therefore, the self-consistency of fj is 
significant for quantitative understanding of electron correlations. 

Within the self-consistent MLA, we have clarified the role of the momentum dependence 
of variational parameters in comparison with the original LA. We demonstrated that the self- 
consistent MLA improves the LA irrespective of the Coulomb interaction energy parameter U 
and electron number n. The correlation energy in the MLA is lower than those of the LA and 
the GA in the weak and intermediate Coulomb interaction regimes. Thus the MLA wavefunc- 
tion should be better than both the LA and the GA in these regimes. The double occupation 
number is suppressed as compared with the LA both in the same interaction regimes. We 
found that the calculated momentum distribution functions show a distinct momentum de- 
pendence. This is qualitatively different from the LA and the GA because both of them lead 
to the momentum-independence of the distributions below and above the Fermi level. We also 
found that the quasiparticle weight in the MLA is lower than that of the LA irrespective of 
U and n, and is close to the result of the NRG. Accordingly the critical Coulomb interaction 
U C 2 of the MLA becomes closer to that obtained in the NRG. 

The variational theories discussed in the present paper construct the correlated ground 
state by applying the two-particle operators to the Hartree-Fock ground state. The wavefunc- 
tions of both the LA and the GA are expressed by the momentum dependent two-particle 
excited states in addition to the Hartree-Fock one. Each amplitude of the excited states 
is momentum independent in these methods. The MLA wave function, on the other hand, 
each amplitude of the two-particle excited states is momentum dependent. By choosing the 
momentum-dependent amplitudes ?7^ fc2 ^ fcl i>est on the basis of the variational principle, we 
improved the LA and the GA in the weak and intermediate Coulomb interaction regimes. Note 
that this procedure does not depend on dimensions of the system within the SSA because of 
the local projection (k[\i) (i\ki) (k' 2 \i} (i\k2) in the operator Oj and thus the improvement re- 
mains unchanged even in infinite dimensions. 

Needless to say, many methods to solve the correlation problems in infinite dimensions have 
been developed. 2 ' 31 ) The NRG is one of the best approaches to calculate the excitations at zero 
temperature as well as related ground-state properties. The accuracy of MLA is on the level 
of the iterative perturbation theory 32 ) at the present stage. Furthermore excitation properties 
cannot directly be calculated by the variational approach. It is also true, however, that the 
high-quality methods such as the NRG 27 ) are not applicable to the realistic systems because of 
their complexity. The present approach is applicable to more complex systems. Moreover the 
wavefunction method arrows us to calculate any static averages as we have demonstrated in 
paper I. 18 ) We therefore believe that further developments of the MLA wavefunction approach 
will provide us with a useful tool for understanding correlated electrons in the realistic systems 
and their physics. Developments of the theory towards the strongly correlated system and its 
extension to the realistic systems are our current problems in progress. 
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Appendix: Laplace transform for the correlation calculations 

The Laplace transform can significantly reduce the number of integrals in the physical 
quantities which appear in our variational theory. It is written as follows: 

= -»/ dte i{z - €4+e3 - e2+ei+ec)t . (A-l) 

z - e 4 + e 3 - e 2 + e 1 + e c J 

Here z = uj + i5, and 5 is an infinitesimal positive number. 

The term (OjOj)o in eq. (39) can be written in the energy as follows: 



4 

EI de r . 

n=l 



n P(<* 

n=l 



/(eit)[l-/(e 2t )]/(^)[l-/(64;)] 



(e 4 - e 3 + e 2 - ei - e c 



U 2 fj 2 lim ! 
*->oJ 



4 

EI de n 

n=l 



n P(€n) 
Ln=l 



/(e lt )[l-/(e 2t )]/(63^)[l-/(e^)] 



(A2) 



(A-3) 



(z - e 4 - e 3 + e 2 - ei - e c ) 2 

Here, e n<T = e n + e a , and e a = eo + [/(nj„ CT )o — /i is the Hartree-Fock level measured from 
the Fermi level p. p(e) in the above expressions denotes the density of states for e/%, i.e., the 
Fourier transform of tij. Now using the relation of Laplace transform (A-l), we obtain 

(dJOi)o = —U 2 fj 2 lim / dtdt'e^ +e ^ t+t '^ / ckie iei (* + *')p(ei)/(ei t ) 

x | d6 2 e- 2 ( t+t ')p(6 2 )[l-/(6 2t )] | d6 3 e^(*+O p(e3)/( ^ ) 

x J de 4 e-^ t+t '^(e 4 )[l-f(e 4l )} 

= -U 2 i? / dtdt'e i€c( - t+t,) a t (-t - t%(t + t') ai (-t - t%(t + i')- 
Jo 



Here 



a a (t) = I dep(e)f(e + e a )e 



-let 



b a (t) = J dep{e)[l- f(e + e a )]e 



ict 



(A-4) 
(A-5) 

(A-6) 
(A7) 



The 4-fold integrals of (OjOi) in eq. (A-2) reduce to the 2- fold integrals in eq. (A-5). 

In the same way, we can perform the Laplace transform of various elements in the physical 
quantities as follows: 

(Hdi) = {6\&)i 



iU 



oo 

2~ / jj. n ie c t 



JO 



dte^a t (-t)a ; (-t)6 t (t)6 4 (t) , 



(A-8) 



(d}H di)o 



-U 2 fj 2 



/ dtdt't 
Jo 



,ie c (t+t') 



[a t {-t - t')b t {t + t') ai (-t - OM* + *0 



-a t (-t - t')b t (t + t')a n {-t - t')b±(t + i') 
+a t (-t - t')b lt {t + 0a4.(-t - t'Mt + i') 
-ai t (-t - t')b t (t + i>+H - t'Mt + 0] , 



(A-9) 
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{d]Oidi) = -U 2 f dtdt / e <ec (* +t, )[a r (-t)6 r (t + Oa4.(- t )^(* + < / ) a r(-* / ) a 4.(-0 

J o 

-a t (-t)6 t (t + f) ai (-t - i>+(i)a t (-t')M*') 
-a t {-t - i>t(*H(-*)M* + t'W^i-t') 
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Here 



a la (t) = J dep(e)f(e + e a )ee- iet , (A-ll) 

6i ff (t) = y dep(e)[(l-f(e + e a )]ee- iet . (A-12) 
The element (34) and (33) for the calculation of the best choose of fj are expressed as 

/•oo 

C = - dtdt'e iec( - t+t,) [a t (-t)b t (t + *>4.(-t)&4.(* + Oat(-*'H("0 

-a t (-t - t')& t (*H(-*)M* + OM*'H(-*') 
-a t (-t)6 t (t + t'HH " 0^(*)«t(-*')64.(* / ) 

+at (_i _ t %(t) ai (-t - t'^b^t'^t')] , (A-13) 



and 



D = i/ dte iec *a t (-t)a ; (-t)6 t (t)6^(t) . (A- 14) 

The correlation contribution to the momentum distribution function (44) is given by 

/•oo 

N(d\h k<T di)o = U 2 f dtdt'e^^a-^-t-t^b-ait + t') 

Jo 

x [f(e k a)ba(t + t>y^ t+t '^ - [1 - f(e ka )]a a (-t - t')^ 1 ^] . (A-15) 

The correlation contribution to the electron number (43) which appears in the calculation of 
the double occupation number is expressed as 

/•oo 

{OjniMo = -U 2 f / dtdt'e ie ^ t+ ^ [a- a (-t - t')b_ a (t + t')a a (-t - t')b a {t)b a {t') 

Jo 

-a- a {-t-t')b_ a {t + t')a a {-t)b a {t + t')a a {t')] . (A-16) 
The element (45) for the calculation of the double occupancy is expressed as 

/•oo 

(d\Oi) + {Oidi) = 2iUfj cfte i£c *a t H)M*H(-W) ■ (A-17) 

J o 
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